UMDGR-94-76 
gr-qc/9312033 



Quantum Gravitational Collapse 
of a Charged Dust Shell 

: Eli HawkinsQ 

^ ■ Department of Physics, University of Maryland, College Park, MD 20742-4111 

o ■ 

Q . Abstract 

I A simple self gravitating system — a thin spherical shell of charged pres- 

sureless matter — is naively quantized as a test case of quantum gravitational 
CS| ' collapse. The model is interpreted in terms of an inner product on the posi- 

[ five energy states. An S'-matrix is constructed describing scattering between 

^ ' negatively and positively infinite radius. 

o : 
^ : 

JTj ■ 1 Introduction 

a^ '■ 

I In the absence of a successful quantum theory of gravity, it is interesting to study 

CT' simplified, naively quantized models of self gravitating systems. Although such a 

model is not backed by any complete theory, it may provide clues as to the properties 
of a complete theory. Such models may also confront some of the poorly understood 
• ^ . issues of quantum gravity in a simplified setting. 

^ I One interesting gravitational process predicted by General Relativity is gravita- 

te I tional collapse beyond an event horizon. A quantitatively accurate model of gravi- 

tational collapse, taking into account the equation of state of the matter as well as 
its rotation and other factors is too complicated to solve in general classically, let 
alone attempt to quantize without a full theory. However, some simple special case 
of collapse might indicate the qualitative behavior of the more general case. In this 
paper, the equation of state is dealt with by assuming the collapsing matter to be 
a pressureless dust. The motion is simplified by assuming spherical symmetry, but 
this still leaves the system with infinite degrees of freedom. If, instead of a solid ball 
of matter, the collapsing object is an infinitesimally thick shell, then it has only one 
degree of freedom, its radius, and it is simple enough to deal with. 

Naive quantization of such a dust shell has been accomplished by Hajicek, Kay, and 
Kuchar [2] . Unfortunately, that model is too simplified to exhibit one of the interesting 
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properties of general gravitational collapse; in general, if the collapsing object has any 
angular momentum or charge it may reexpand, possibly into an asymptotically flat 
region distinct from the one in which it originally collapsed. Allowing the shell to 
have angular momentum would break its spherical symmetry and complicate matters 
greatly. A charged shell, however, would exhibit some of the same phenomena, but 
could still be spherically symmetric and have only one dynamical degree of freedom. 

Unfortunately, a major qualitative difference remains between the charged and 
the rotating cases. In the Kerr space-time there is an accessible, asymptotically flat 
region where r < 0; in the Reissner-Nordstrom metric such a region exists, but it 
is only connected to r > at a single point of space, and thus hardly seems like 
a physical part of the space-time. It seems plausible nonetheless to assume this 
spherically symmetric model to be continuous with the more general case of nonzero 
angular momentum and therefore to take the negative radius region seriously. 

In this paper I naively quantize the collapse of a charged dust shell using an 
equation of motion that is a trivial generalization of the Wheeler-DeWitt equation 
derived in [2]. However, instead of applying a boundary condition at the origin and 
only allowing positive radii, I apply a matching condition at the origin and allow the 
shell to attain negative radii. 

One issue which this model cannot be expected to resolve is whether singularities 
occur in quantum gravity. If there is a singularity avoidance mechanism, then it may 
well involve nonradial degrees of freedom — even in a spherically symmetric situation; 
but all nonradial degrees of freedom are suppressed in this model. 



2 Equations of Motion 

Consider a spherically symmetric, infinitesimally thick, dust shell of rest mass M and 
charge Q. The interior of the shell is a piece of Minkowski space-time; the exterior 
is a piece of Reissner-Nordstrom space-time. The mass appearing in the Reissner- 
Nordstrom metric is the total energy (E) of the system as determined by a distant 
observer, the charge is that of the shell. Note that Schwarzschild coordinates on the 
interior and exterior do not join smoothly across the shell {grr is discontinuous), but 
this does not cause any problem in deriving the equation of motion. 

The classical equation of motion was derived by Kuchaf [1] from the relationship 
between the surface energy tensor of the shell and the difference between the extrinsic 
curvatures of its interior and exterior surfaces. The shell's motion is described by the 
flrst integral of motion 

where R is the radius of the shell, ^ = T is the well deflned Minkowskian time 
coordinate on the interior, and k = M'^ — Q'^. Units are such that h = c = G = 1. This 
equation of motion is formally identical to that of a relativistic point particle moving 
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radially in a Coulomb potential in flat space-time; this analogy is useful because it 
suggests a ready made interpretation for the quantum model. 

By assuming the externally measured energy E to be the energy conjugate to T, 
numerically equal to the Hamiltonian that generates motion in T, the momentum 
conjugate to R is found to be 

r.R MR 

= _ (2) 

and from this a super-Hamiltonian 

is found which is constrained to vanish and which generates correct equations of 
motion. See [2]. 

Promoting the momenta in h to operators makes the constraint h = into a 
Wheeler-DeWitt equation, namely 

Consistent with the above analogy, (^) is equivalent to the Klein-Gordon equation 
for a scalar particle in a Coulomb potential under the replacement = i?$ (equivalent 
to a change of measure). Equation (Hf) is the basis of this naive quantization. 



3 Probability Interpretation 

This quantum model is most appropriate to an observer located in the interior of the 
shell. T and R are well defined in terms of an observer at r = who measures the 
shell's motion by timing the delays of light pulses reflected off the inside of the shell. 
This model may not, however, be able to describe observations made by an observer 
outside the shell. An important effect detectible only by an external observer would 
be the shell passing beyond an event horizon. In this model there is no qualitative 
distinction between the motion of the shell outside and inside an event horizon, as 
there would be in a model based upon an exterior observer who only sees R > 2E. 
This is emphasized by the invariance of both the classical and quantum equations of 
motion under the rescaling: 

(T, R, M, E, Q) — > {aT, aR, aM, aE, ^ - (1 - a'^)M^^ (5) 

for any a > 0. This equates points inside the event horizon in one case with points 
outside the event horizon in another. 
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Following the scalar particle analogy Hajicek et al [2] construct "charge" (really 
number), and energy currents which are generalized as locally conserved bilinear 
currents: 



^^(^1,^2) 
^''(^1,^2) 



2 

i 

2 
1 
2 

-1 



^1^2 -^1^2] +^^1^2 



^^'1^2 + '^''1^2 



(6) 
(7) 



2R 



(9) 



where ' denotes 

Assuming J(\1/,\I^) to be a probability current gives a tentative interpretation of 
the wavefunction for > 0. In the defining experiment for this model an observer at 
r = sends a light pulse at time U which reflects off the shell at (T, R) and is seen 
again at time V. T and R are given by 



T = \{V + U) 
R = \{V-U). 



(10) 



In the quantum model, if the light pulse is sent at the observer's time U then the 
probability that it will return at time V is ^J^dV where = — J^, evaluated 
at {T,R) given by (|10]). (I'm assuming an appropriately normalized wavefunction.) 

This interpretation fails, however, because the probability given is not necessarily 
positive. Following the scalar particle analogy further would suggest that the nega- 
tive densities might represent pair creation, but this model is too crude to make a 
meaningful prediction of that sort. It would be like modeling a car crashing into a 
wall as a potential barrier, and concluding that automobiles would be pair produced. 



3.1 Boundary Condition 

If only i? > is allowed then J and E will only be globally conserved if a boundary 
condition is imposed. Namely, 

J^(^l,^2)=^''(^l,^2)=0 (11) 

at R = for all T. 

This boundary condition rules out ab initio a truly irreversible collapse. When 
the shell reaches radius it is required to pass through itself and reexpand. However, 
the collapse may appear irreversible; if the shell reaches = it is within an event 
horizon and an external observer would never see it reexpand. 
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Negative radius might seem unphysical, but as I have mentioned there are plausi- 
ble reasons for considering it. Among the stable, approximately spherically symmetric 
vacuum space-times (i.e: black holes) the Reissner-Nordstrom space-time is excep- 
tional. For any such space-time with nonzero angular momentum the causal structure 
is that of the Kerr space- time, and there is an accessible region where r < 0. In 
the quantized case, classical states are smeared and the Reissner-Nordstrom causal 
structure might well give in to its generic neighbor. 

In view of the interpretation in terms of an observer at r = 0, allowing R < 
seems especially meaningless — an observer could not possibly continue to exist inside 
the shell after it has completely collapsed. However, if the shell rccxpands then 
observations could be recommenced and their correlation with previous observations 
would depend upon the R < behavior. 



3.2 Revised Interpretation 

When all values — oo < R < +oo are allowed, J gives a conserved inner product equal 
to the flux of J across any Cauchy surface; that is: 

/ + 00 
J^{^i,^2)dR. (12) 
-oo 

Because J^(^, ^) is not always positive, it is not clear that this inner product is 
positive definite, as is necessary for normalizable wavef unctions. In fact it is positive 
definite if and only if the wavcfunctions arc restricted to the subspace of positive 
energy solutions (see appendix). This restriction seems acceptable and excuses the 
failure of the initial interpretation. The local currents were not meaningful because 
the position operators i? or y do not leave the positive energy subspace invariant and 
are thus not good observables. 



4 Classical Motion 

First, consider the motion through i? > 0. In [3] Boulware analysed the motion of a 
charged shell through the possible regions of the extended Reissner-Nordstrom space- 
time that its exterior is a piece of. The Schwarzschild coordinates do not completely 
specify the position of the shell. The key tool which helps to distinguish between 
regions with identical coordinates is the r component of the unit outward normal 
to the outer surface of the shell, which points from the interior (Minkowski) to the 
exterior (Reissner-Nordstrom) region. This is given simply by: 

, E M' + Q^ 

n — . 13) 

M 2MR ^ ' 

When > the exterior of the shell is the r > R region of Reissner-Nordstrom 
space-time. When rf < the "exterior" of the shell is the r < R region of Reissner- 
Nordstrom. 
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For sufficiently small R > 0, n'^ is necessarily negative. This means that the shell 
cannot collapse to form a singularity of charge Q. If it collapses to i? = it in fact 
collides with a singularity on its exterior with charge —Q. 

Setting i? = in and solving for R gives a turning radius 

R - ^ _ M^-Q^ 
+ 2{M-E) 2{M-E)' ^ ' 

Consideration of rf' and R+ gives the principal distinct classical trajectories 
through R> (along with trivial cases such as /c = 0). 

1: For E > M and R+ > (hence \Q\ > M), the shell starts from R = oo, falls 
to a finite radius and bounces back to oo. A further distinction is 

a) \Q\ > E and no event horizon occurs, or 

b) \Q\ < E, an event horizon is formed and the shell reexpands in a different 
asymptotically flat region than it collapsed from. 

2: For ^ > M and R+ < (hence \Q\ < M), The shell either falls from R = oo 
to 0, or the time reverse, exploding from to oo. 

3: For E < M, the shell cannot reach R = oo. Either, 

a) E < ^ and the shell never emerges into an asymptotically fiat region, 

or 

h)E> the shell emerges from an event horizon, expands to a maximum 

radius of then recoUapses, hitting singularities in its past and future. 

Because of the rescaling symmetry (|^), there is no qualitative distinction between 
the solutions for the a and h cases. An interior observer can only guess if the shell 
has passed within an event horizon. 

4.1 Extension to < 

Transforming r —r in the Reissner-Nordstrom metric shows that the r < domain 
of the complete Reissner-Nordstrom space-time is an asymptotically fiat region con- 
taining a singularity with mass — i?, and charge either ±(5; by following the electric 
field lines through r = the charge is seen to be —Q. From (|1^) rf must be positive 
when i? < 0; therefore, if the shell enters i? < 0, this singularity is on the exterior of 
the shell within a compact region. The r component of the normal to the inner surface 
of the shell also changes sign; the interior thus continues to be a ball of Minlcowslci 
space. 

If the classical solution is followed past i? = 0, the radical in (||) must change sign; 
this analytically continued equation of motion is presumably the only one continuous 
with that for nonzero angular momentum. There is now a second turning point at 

R -fc _ M^-Q^ 

2{M + E) 2{M + E)' ^ ' 
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For case 2, the shell passes in a finite proper time to R = R- < and then bounces 
back out to R = +00. For case 3, The shell oscillates between i?+ > and i?_ < 
with a finite period. So, classically, if the shell is allowed to pass into < 0, it will 
always return to -R > and potentially be observed again, except in the special case 
= in which the equation of motion reduces to i? = Const. 



5 Solution of the Wave Equation 

Now find a complete set of solutions to the Wheeler-DeWitt equation (|) by consid- 
ering definite energy states {P^ eigenstates) \E'(T, R) = e~''^'^'i/j{R). With this ansatz, 
(I) becomes 



= ^" + 



E+ — ] -M' 
. 2R 



,^n\ Ek k"^ 



Now assume E > M. Defining 



X 



Ve^ - M2 ; x = xR 



ip. (16) 



(17) 



becomes 







+ 



Ek fc2 
XX Ax"^ 



The characteristic equation about the regular singularity at x 



is 



= A(A-1) + 



4 



:i8) 



(19) 



with roots X± = where a = \/l — k'^. As noted in [2], this model breaks down 
for I A; I > 1, corresponding to an overcritical Coulomb potential which can only be 
described by a second quantized theory. It is far from clear what an appropriate 
analogue of second quantization would be here, and it also seems implausible that an 
analysis of higher order effects could be meaningful in a model with virtually all the 
system's degrees of freedom suppressed. I therefore only consider \k\ < 1. 

First, only consider the solution of (|18D in x > 0; there must exist independent 
solutions such that 

^ = (1 + 0{x)) . (20) 
Because the coefficients in (|T^) are real, the solutions (|20|) are real for all x > 0. 



Adopting one more new constant f] = ir-, (pTSf) may be written 



d^ij 

d{2ix) 



+ 



2X 

1 

'4 



2ix 



+ 



4 4" 



{2ix) 



(21) 
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which is just Whittaker's equation (see [4] p505). Define two functions, equal to the 
solutions (pop: 

V?±(a, p;x>0) = e^"x^±M (A± + i(3, 2A±, 2ix) (22) 

where M is the regular Kummer confluent hypergeometric function. These solutions 
are indeed real (see [4] eq. 13.1.27). 

M is a single valued function, but because of the x^^ factor ip± are generally 
multi-valued functions in the complex plane with a branch point at x = 0. In order 
to analyse the R < behavior it must be determined which branch the physical —x 
axis lies on. (This analysis also applies when E < M.) As a criterion, note that (|^) is 
invariant under the time reversal transformation \l/(r, R) —>■ ^{—T, R) or equivalently 
ip{R) i^{R), E ^ E. Because of this we will demand that the complex conjugate 
of any physical solution of (^) must also be a physical solution. This means that Lp± 
must be real for x < as well as x > 0, otherwise ip± and (p± would constitute more 
than 2 independent solutions to (|16D . 

A second criterion is provided by the current form J^. For two deflnite energy 
states of the same energy, conservation implies that must be constant with respect 
to R. So 

(^+(a, /?; x)ip'_{a, (3; x) — f'^{a, /?; x)ip^{a, /?; x) = Const. (x). (23) 

This implies that one of the (p± is approximately even, and the other approximately 
odd near x = 0; it doesn't matter which. Arbitrarily choose as approximately 
evenQ. 

Now note that (pTf ) is invariant under {x,(3) (— x, — /3) so 

Lp±{a,p;x) = ±Lp±{a,-p;-x). (24) 



6 Transmission and Reflection coefficients 

In a definite energy state for R ±oo, J"^(\I^,\I^) becomes simply E ■ l^^f; the 
interpretation simplifies. For a fixed energy, the wavefunction near infinity may be 
interpreted as a relative probability amplitude, and an S'-matrix may be constructed. 
Time reversal invariance and unitarity restrict the S- matrix to the form: 




S={ 7^ 7- ) (25) 

where 



5"++ = e*'' cos 7 
S^^ = ie*'^sin7 



iS.;^., (26) 



^Strictly speaking, this sign combination (or its alternative) cannot be attained as a branch of 
the multi- valued solution for certain rational values of a, but the complement of that set of values 
is dense, and this is physics. 
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If (^) is viewed as a time independent Schrodinger equation, the effective poten- 
tial is long ranged, and so the wavefunctions cannot asymptote to definite momentum 
waves with stable phase. As a result the phase 6 is not well defined; only the reflec- 
tion and transmission coefficients cos^ 7 and sin^ 7 are meaningful. The reflection 
coefficient is interpreted as the probability that a shell collapsing from +00 will re- 
turn to +00; assuming this model to have any validity, that probability should be 
independent of how the position of the shell is measured. 

These reflection and transmission coefficients for the scattering of wave packets 
between -|-oo and —00 are determined from the asymptotic behavior of the two basis 
solutions ijj{R) = ip±{a, l3]xR)- 

Equation 13.5.1 of [4] gives the asymptotic behavior of Lp±: 



'fi±{a, (3;x^ +00) = (±a)! 
and by (H), 

ip±{a, [3]X ^ —00) = ±(±a;) 



i \ 
2 



) r(A±-^/?) 



-|- c.c. 



(l + 0(x-i)) (27) 



-2x) 



-if) 



r(A± + ip) 



+ c.c. 



(1 + (x-i 



f28) 



In lieu of definite momentum waves, the appropriate ingoing and outgoing asymp- 
totic waves are e^""^ where 



(P = xR + sgn{R)(3\n\2xR\ 
In terms of (p, the asymptotic forms of the wavefunctions are 



ij{R +00) = (±a) 

i;{R-^ -00) = ±(±a): 
which are of the form 



--0 



r(A± - ti3) 



e"*^ + c.c. 



(1 + [r-^ 



-0 



r(A± + ti3) 



e''^ + c.c. 



1 + i? 



+00) 
—00) 



(29) 

(30) 
(31) 

(32) 



By applying the definition of the ^-matrix: 



CoutH 
Cout - 



s 



Cin+ 
Cin — 



one obtains, 



±S. 



^4 



r {x± - 1(3) 



(33) 



(34) 
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Combining the + and — equations and solving gives 



6 = arg [r (A± - i(3) T (A± + i(3)] (35) 



TT 



cos^7 = sin^ sech^irf]. (36) 

Classically, the shell cannot reach R = —oo, unless = in which case it must. 
Correspondingly, the transmission coefficient, cos^7, is necessarily < 1 unless k = 0, 
in which case it is necessarily = 1. The transmission coefficient goes to when 
k ^ ±1 or E ^ M. 

There are poles in the transmission coefficient when coshyr/? = 0; accepting only 
normalizable (positive energy) solutions, this implies the existence of bound states 
with energies 

E„ ^ /^<^"-^> (37) 



{2n - 1)2 + P 

with n=l,2,3,.... This is a different energy spectrum than results if the boundary 
condition ( pJ]) is imposed. For A; < 0, the bound states correspond to states with no 
asymptotically flat region. For k > the bound states all correspond to classical states 
with maximum radius (-R+) outside the event horizon (classical case 36). Because this 
model is unable to deal with an observer outside the shell, it cannot decide questions 
of when the shell does or does not pass beyond an event horizon any better than this. 

The bound state solutions can be found explicitly in the same way as the scattering 
solutions, but in contrast to the bound states found in [2] with the boundary condition 
(P!T|) imposed, they are not algebraically special. 



7 Conclusions 

By necessity, this Quantum model is based upon several tentative assumptions. The 
standard Schwarzschild time t is the proper time of a distant, stationary observer; 
such an observer will measure the total mass of the system to be E; therefore E is 
the energy conjugate to t. The ffist assumption made in this model was that E is 
conjugate to T. This works if there is a time slicing such that T and t differ only 
by a constant; that may be accomplished locally but not necessarily globally in this 
system. This issue will be further addressed for the uncharged case in [5]. If the 
correct energy for T is not E, then it must be some constant of motion invariant 
under time translations (as commutes with itself); this restricts it to be some 
homogeneous function of E, M, and Q. 

The second assumption is that a Wheeler-DeWitt equation based on the super- 
Hamiltonian is the correct method of quantization; this is probably not completely 
true and is the naive aspect of this model. 
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The third assumption is that, based on the equivalence of the equations of motion 
with those of a scalar particle, the inner product (|12D gives the correct probability 
interpretation of the wavefunctions. 

The restriction of my interpretation to measurements made by an observer on the 
interior is necessary. One might attempt to apply this model to observations made 
from oo; if a distant observer reflects light off of the shell (to be precise a pulse is 
sent from some X~ and received at some X+, possibly in a different asymptotically 
flat region), then the appropriate coordinates are the Kruskal coordinates of the 
shell in the exterior region. Unfortunately, the transformation to these coordinates 
from T and R depends upon the whole trajectory of the shell; there is no obvious 
quantum equivalent. A model analogous to this one, but based upon an outside 
observer appears intractable. The first step to constructing it would be to write E 
as a function of R and but this is the solution of a quartic equation and is rather 
unwieldy. 

This restriction to interior observers means that this model cannot directly deal 
with the interesting questions of when (or rather with what amplitude) the shell will 
pass beyond an event horizon and remerge in a different asymptotically fiat region. 
This can only be addressed crudely by comparing quantum states with classical tra- 
jectories. A typical wavepacket will be a combination of definite energy states whose 
classical equivalents do and do not contain event horizons. This combination can be 
crudely interpreted to give the probability of an event horizon occurring. 

This model has in some sense addressed the issue of whether singularities occur. 
When i? < 0, the corresponding classical space-time will include a singularity at 
r = 0. However, the dynamics of this model are only mildly singular at i? = 0, 
mildly enough that the solution could be simply continued across that singularity. 
This might be taken as a hopeful sign that the dynamics of a full quantum theory of 
gravitation will be singularity free. 



8 Appendix 

Proposition: The restriction to positive energies makes the inner product ^T^ ) positive 
definite. 

To see this consider a state \E' = e~^^^il){R) with definite energy E' > (it is 
straightforward to show that different energy states are orthogonal). We have 



For scattering states this must be positive as |-R| ~ oo contributions dominate. For 
bound states ip^R ±oo) — > 0, hence for a > 0, 







{^{aR))'^{R) - ^{aR)%lj\R) 



+ 00 

— oo 
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/- 



— oo 



(V^(ai?))Xi?) 



ip{aR)ilj"{R) dR 




p+oo / \ _ 

{I -a) [1 + a]^ + - - (1 + a)M2 tP{aR)^{R) dR, 



J —oo \ R 



so for a 7^ 1 



;>+oo / \ _ 

= / E [[l + a]E+ - ] - {1 + a)M^ i;{aR)ip{R)dR. 



J —oo \ R 



(39) 



This integrand diverges at i? — > as i?^" (recall |a| < 1) and approaches exponen- 
tially as — > ±oo. As a result, this improper integral may be uniformly approximated 
to arbitrary accuracy by a proper integral and is therefore a continuous function of a 
and so 



which is manifestly positive. Similarly this is clearly negative ior E < 0; hence 
restriction to positive energy is both necessary and sufficient for the inner product to 
be positive definite. 
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